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Abstract 

The incorporation of Phosphorus into (Ga,Mn)As epilayers allows for the tuning of the magnetic 

H ! easy axis from in-plane to perpendicular-to-plane without the need for a (Ga,In)As template. For 

'^ ' perpendicular easy axis, using magneto-optical imaging a self-organized pattern of up- and down- 

Sin ; 
O ■ magnetized domains is observed for the first time in a diluted magnetic semiconductor. Combining 

Kerr microscopy, magnetometry and ferromagnetic resonance spectroscopy, the exchange constant 

K^ I and the domain wall width parameter are obtained as a function of temperature. The former quan- 

0^ ■ tifies the effective Mn-Mn ferromagnetic interaction. The latter is a key parameter for domain wall 

dynamics. The comparison with results obtained for (Ga,Mn)As/(Ga,In)As reveals the improved 

in 

^D ' quality of the (Ga,Mn)Asi_„Py layers regarding domain wall pinning, an increase of the domain 

o ■ 

wall width parameter and of the effective Mn-Mn spin coupling. However, at constant Mn doping, 
no significant increase of this coupling is found with increasing P concentration in the investigated 
range. 



X 



PACS numbers: 75.50.Pp, 75.60. Ch, 75.70.Ak 



Diluted ferromagnetic (FM) (III,Mn)V semiconductors are the subject of intensive re- 
search in view of potential applications in the field of spintronicsi'^. Several proposals 
for magnetization manipulation using either light-induced magnetization switching between 
different magnetic easy axes^^^ or magnetic domain wall (DW) motion^"^ have emerged. 
(Ga,Mn)As layers and microtracks with perpendicular magnetic anisotropy are well suited 
for the investigation of DW propagation. In such layers the intrinsic flow regime for field- 
driven DW propagation was demonstrated for the first time in a diluted FM semiconductor-. 

One of the properties of this versatile FM semiconductor is the strong dependence of 
its magnetic anisotropy on epitaxial strain^'^.-ii. This dependence results from the carrier- 
mediated origin of the ferromagnetism. It arises from the competition between the Mn- 
induced giant Zeeman splitting and the strain-induced splitting of the valence bands. Re- 
cently, it was shown that introducing a few percent of Phosphorus in substitution of Arsenic 
in (Ga,Mn)As layers grown on GaAs provides an efficient tuning of the epitaxial strain, 
from compressive to tensile, due to the smaller ionic radius of Phosphorusi^ii^. As expected 
theoretically, the magnetic easy axis switches from in-plane to out-of-plane orientation upon 
increasing P concentration. This eliminates the need for a highly mismatched (Ga,In)As 
buffer layer for realizing out-of-plane anisotropy. This buffer layer was shown to introduce 
a cross-hatch pattern and emerging dislocations which are detrimental to the propagation 
of magnetic DWs^'ii. (Ga,Mn)Asi_j^Pj^ films grown directly on a GaAs buffer are expected 
to contain much less defects. In addition to this structural advantage, (Ga,Mn)Asi_j,Py 
films are appealing for two aspects-^. Firstly, it is expected that the smaller lattice con- 
stant of (Ga,Mn)Asi„j^Pj, will hinder Mn incorporation in interstitial sites, where Mn ions 
act as charge and magnetic moment compensating defects, thereby favoring higher carrier 
density than in (Ga,Mn)As. Secondly, it was predicted that the decrease of the lattice con- 
stant should also lead to an increase of the Mn-hole exchange integral Jpd and hence of the 
Curie temperature. In turn, this should lead to an increase of the exchange constant A 
describing the effective exchange interaction between the Mn spins in the framework of the 
micromagnetic theory. The determination of A yields the DW width parameter defined as 



A = ^J AjKu-, where K^ is the uniaxial anisotropy constant. A is a key parameter for DW 
dynamics^-. The Walker velocity, i.e. the theoretical critical velocity at the crossover from 
the steady to the precessional flow regimes in the one-dimensional model, depends on A as 
'^w = 7A/io^s/2, where 7 is the electron gyromagnetic ratio and Mg the magnetization. 



Furthermore, in both the steady and precessional regimes, the DW mobihty is proportional 
to A. 

In FM layers with perpendicular easy axis one of the methods for the determination of 
the exchange constant A is based on domain theory, i.e. micromagnetic theory applied 
to self-organized magnetic domainsi^"— . Self-organization of up- and down- magnetized 
domains in a periodic array results from the competition between the DW energy and the 
magnetic energy arising from long-range interaction between domains. Up to now, although 
an estimation of A could be obtained from the domain structure and hysteresis cycle in 
(Ga,Mn)A3i^, a long-range periodic arrangement of magnetic domains in (III,Mn)V FM 
semiconductors had not yet been obtained. 

In this paper we report the observation of self-organized magnetic domains in (Ga,Mn) Asi^^^Pj^ 
FM alloys with perpendicular magnetic anisotropy. We determine the exchange constant A 
and the DW width parameter A from the domain period. An effective exchange coupling 
constant JmuMu between Manganese spins is then obtained as a function of Phosphorus 
concentration. 



The samples were 
can be found in Ref. 



gown by molecular beam epitaxy on GaAs (001) substrates. Details 
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131]. For (Ga,Mn)Asi_,yPj^ samples, the Mn (P) concentration was 
estimated from a reference (Ga,Mn)As ((Ga,As)P) sample, grown under similar conditions, 
in particular a similar substrate temperature. However the determination of Mn and P con- 
centrations is made difficult by the presence of an unknown concentration of interstitial Mn 
and As in antisites. In sample D, with no Phosphorus, the (Ga,Mn)As layer was grown in the 
traditional manner on a Gai-^In^As buffer layer (z=0.098) in order to achieve perpendicular 
anisotropy^. To improve M^, the Curie temperature Tq and the carrier concentration the 
samples were annealed under N2 atmosphere at 250° C for 1 hour. For all samples the layer 
thickness is (i = 50 nm. 

The magnetic domain structure was investigated using polar magneto-optical Kerr 
(MOKE) microscopy at variable temperatureH. The lateral spatial resolution is 0.9 yum. 
The light wavelength is 600 (670) nm for the (Ga,Mn)AsP ((Ga,Mn)As) samples. The hys- 
teresis cycle was obtained from the average intensity of MOKE images as a function of the 
applied field. The temperature dependence of the saturation magnetization Ms{T) and the 
Curie temperature Tc were obtained from magnetometry using a superconductor quantum 
interference device (SQUID). The anisotropy constants were determined from ferromagnetic 



resonance (FMR) spectroscopy. The spectra were analysed using the Smit-Beljers equation 

and the minimization of the free energy for different ahgnments of the apphed magnetic 
fieldi^'20. 

Table [I] summarizes the characteristics of the six samples which were selected in order to 
compare layers with the same Mn concentration and different P concentrations (A series, C 
and D), or the same P concentration and different Mn concentrations (A2 and B, A3 and 
C). Most of the samples show a Tc around 110 K. Tc is not improved significantly in the 
(Ga,Mn)Asi_yPy samples as compared to sample D (no Phosphorus). The concentration of 
Mn ions participating in the ferromagnetism [Mue//] is determined from the saturation mag- 
netization Ms at T = 4 K as [Mneff]=Ms/{NQgfiBS), where A^o is the density of cation sites 
in the lattice, fiB the Bohr magneton and S = 5/2. For all samples the ratio [Mne//]/[Mntoiai] 
is ~0.5, revealing the fact that the incorporation of Mn in interstitial sites is not modified 
by alloying with Phosphorus. 

The method employed to determine the exchange constant and the DW width parameter 
is based on domain theory for layers with perpendicular easy axis^^. Self-organization of 
up- and down-magnetized domains in a periodic stripe array minimizes the total energy 
consisting of the Zeeman energy, the surface energy of the DWs and the stray field energy. 
The stripe period p (in units of d) in the demagnetized state (zero applied field and zero 
average magnetization) satisfies the following equation^^: 



p oo Sin'' (f ) X ^-pQ - 2mT{l + Q) + pQ (Cosh (x) + ^l + ^Sinh (x) 



^^^ n=i n3 (Jl^Cosh (x) + Sink (x) 



where x is equal to i2m:^/l + l/Q\ /p. Q = Ku/{fioMg /2) is the anisotropy quality factor 
shown in Table [B Q is smaller in the P incorporated samples than in sample D due to the 
smaller K^ and larger Mg. Ac is Thiele's length defined as a/fioM^d with a the specific 
DW energy. Knowing p and Q, one determines Ac and hence a. Typically, Ac is found to 



be of the order of 1. For a Bloch wall a is equal to Ay/AKu- Micromagnetic simulations 
using OOMMF package^ show that in our samples with d=50 nm DWs are not Bloch walls 
but twisted walls with Neel caps. Using the variational calculations of Ref. [22j for the 
DW specific energy of twisted walls we obtain the functional dependence of a on Q and A 
namely a = 4:f{Q, A)AKu with f{Q, A)< 1. For each value of a we obtain the DW width 



parameter A and hence the exchange constant A. 

Figure [T](a) shows a typical MOKE image of the self-organized domains with up- and 
down- magnetization (white and black) after a demagnetization process in an alternating 
magnetic field of decreasing amplitude. The remnant field is < 1 G. This image is obtained 
on sample A3 at T = 80 K. Figure [TJ^b) shows the domain period obtained from the MOKE 
images as a function of the temperature T for this sample. The period is found in the 
range 2-3 /xm for 4 K< T <105 K. The period for 4 K< T <40 K and T=105 K is 
determined from MOKE images with a weak optical contrast. We will discuss below the 
reliability of the data extracted from these values. Following the procedure described above, 
Thiele's length Ac, the DW specific energy a, the DW width parameter A, and eventually the 
exchange constant A are obtained. A and A are shown in Fig. [2] for sample A3 (filled square 
symbols) as a function of temperature. A is found in the range 7-9 nm and A in the range 
0.1 — 0.4 pj m~^. The decrease of A with temperature mainly reflects the decrease of the 
magnetization {A oc X'^M'^/Ku)- In order to assess the reliability of the A and A values at 
low and high temperature, we determine the upper and lower boundaries for these quantities 
using the method previously developed for (Ga,Mn)As/GaInA3iS. The upper boundary for 
Ac and consequently for A and A is determined from the width of the stripe domains of the 
minority phase observed close to the saturation field. For sample A3 this upper boundary 
coincides at all temperatures with the values obtained from the domain period. The lower 
boundary for Ac (and hence for A and A) is obtained from the comparison of the normalized 
experimental hysteresis cycle with the calculated magnetization curve m{h), where m is 
the spatially averaged reduced magnetization (M) /Mg and h the reduced applied magnetic 
field H/Mg. The m{h) curve is obtained from the minimization of the free energy of the 
periodic stripe array. Ac is the only free parameter. An example of experimental hysteresis 
cycle and calculated m{h) curve is given in Fig. |3l The m{h) curve which is tangent to 
the experimental hysteresis cycle provides the lower boundary for Ac and hence for A and 
A. The values obtained by this method are shown by open symbols in Fig. [2l At high 
temperature (90-100 K) the lower boundary values coincide with the values determined 
from the domain period. To summarize, between 4 K and 40 K the DW width parameter 
and exchange constant are comprised between the values shown by the closed and open 
symbols of Fig. 2, respectively. Above 40 K A and A are reliably obtained from the domain 
period (filled symbols). Similar results are obtained for samples Al and A2 where the self- 



organized domain pattern can be achieved. For samples B and C DW pinning by defects 
prevents the self-organization of domains in regular patterns. Therefore only the lower and 
upper boundaries of A and A can be obtained, using the method recalled above. For sample 
C the upper boundary values could not be obtained below 30 K since the stripe domains 
were not observed. For sample D ((Ga,Mn)As on a (Ga,In)As template) the DW width 
and exchange constant were obtained from the analysis of the field-driven DW dynamics^ in 
excellent agreement with the estimation from the domain structure and hysteresis cycled. 

A direct comparison of the exchange constant A between the A, B, C, and D samples 
is not relevant because of their different effective Mn concentrations. Therefore, using Kit- 
tel's relation between the first-neighbor spin-spin interaction constant and the exchange 
constant A of the micromagnetic theory^, we obtain an effective interaction constant be- 
tween Manganese spins JunMn independent of the Mn concentration as JunMn = Aa/2S^ 
where a = {2/Nq [Mneff]Y^^ is the spin lattice parameter—. Figure H] represents JmuMu as 
a function of P concentration. The JmuMu values are shown for T /Tc ~ 0.4, which is the 
lowest temperature for which data are available for all samples. JmuMu does not scale like 
Tci which points to the importance of distant Mn interaction via the hole gas^^. JmuMu is 
larger for the P incorporated samples than for the (Ga,Mn)As sample by a factor up to 8. 
However, among the A series there is no clear tendency toward an increase of JmuMu with 
P concentration, at least within the narrow range investigated here. These samples with 
the largest Mn concentration (10.4%) also have the largest JmuMu, which may result from 
a larger hole concentration. These results suggest a possible effect of Phosphorus incorpo- 
ration on the Mn-hole exchange integral Jpd- A quantitative estimation of Jpd is however 
not possible without the precise determination of the carrier concentration and a theoretical 
description of the hole density of states in the (Ga,Mn)Asi_yPy alloys, which is beyond the 
scope of this paper. 

Figure [5t^a) shows the DW width parameter A as a function of P concentration. A is 
larger in P incorporated samples. This suggests that the DW mobility (field derivative of 
the DW velocity) in the field-driven DW propagation regimes, which scales like A, should be 
larger in the P incorporated samples. We then obtain the Walker velocity vw = 7A;UoMs/2. 
We predict an increase of Walker velocity for the P incorporated samples, as shown in 
Fig. E^b). DW propagation experiments are needed to confirm these predictions. 

As a conclusion, the novel ferromagnetic semiconductor (Ga,Mn)AsP is a very promising 



material for achieving FM layers with perpendicular easy axis. We have shown that the 
determination of the exchange constant and DW width parameter from domain theory is 
made more reliable than in (Ga,Mn)As owing to the self-organization of magnetic domains. 
This is made possible by the low density of DW pinning centers. 

From the exchange constant A an effective interaction constant between Mn spins 
JmuMu has been obtained. (Ga,Mn)Asi_j,Pj^ samples have a larger JMnMn than the 
(Ga,Mn)As/GaInAs sample, however no significant enhancement of JmuMu as a function of 
P concentration is found in the range investigated here. The Curie temperature and the 
ratio [Mneff\/[Mntotai\ do not exhibit any remarquable improvement. 

The DW width parameter A as well as the predicted DW Walker velocity show a signif- 
icant increase in the P incorporated samples. This should be of great importance for the 
study of field-driven and current-driven DW propagation. 
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sample 


Al 


A2 


A3 


B 


C 


D 


[Mnt„t,,](%) 


10.4 


10.4 


10.4 


8 


7 


7 


[P](%) 


11.3 


8.8 


7 


8.5 


7 





Im (ppm) 


-6650 


-5120 


-4260 


-6870 


-3580 


-10700 


ezz(%) 


-0.3 


-0.2 


-0.2 


-0.3 


-0.2 


-0.5 


Tc(K) 


119 


113 


139 


110 


80 


125 


Ms (kA m-i) 


53.8 


50.8 


53.5 


46.2 


39.2 


38.5 


Q 


4.35 


3.3 


3 


4.8 


3.25 


8.5 
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FIG. 1. (a) Self-organized pattern in the demagnetized state for sample A3 at T = 80 K. (b) 
Domain period as a function of temperature for sample A3. 
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FIG. 2. Exchange constant A (a) and DW width parameter A (b) for sample A3. The filled and 
open symbols represent the values determined from the period of self-organized magnetic domains 
and from the hysteresis cycle, respectively. 
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FIG. 3. Experimental hysteresis cycle (symbols) for sample A3 (10.4 % Mn, 7 % P) at T =100 K 
and calculated ?7i(/i) curve for Ac = 0.95 (solid line). 
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FIG. 4. (Color online) Effective exchange energy between the Mn spins JMnMn as a function of P 
concentration at T/Tc ~ 0.4. For samples B and C the upper and lower boundaries for JMnMn 
are shown by filled and open symbols, respectively. 
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FIG. 5. (Color online) DW width parameter A (a) and Walker velocity for DW propagation (b) as 
a function of P concentration at T/Tc ~ 0.4. For samples B and C the upper and lower boundaries 
for these values are shown by filled and open symbols, respectively. 
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